function model=update_zeta_1(model)


temp      = (sum(exp(model.alpha_1 + (1/2)*model.sigma_1.^2), 2))';
ind       = find(abs(temp)<model.MINVALUE);
temp(ind) = model.MINVALUE;

model.zeta_1 = temp;

end
